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1 .  INTRODUCTION 

This  paper  investigates  alternatives  to  minimum 
variance  unbiased  (MVU)  estimators  in  noncentral  (NC) 
X2  and  P  distributions.  Motivation  is  provided  by 
noting  that  if  Z  ^  X2(Ps*),  then  i)  (Z-p)/2,  the  MVU 
for  g(A)  =  A  is  inadmissible  under  squared  error  loss 
(SEL),  see,  e.g.,  Perlman  and  Rasmussen  (1975)*  ii) 
Z"1,  the  MVU  for  g(A)  =  E(Z-1)  is  inadmissible  under 
any  "bowl  shaped"  loss  since  S(Z)  =  min(Z  1,(p-2)_1) 
dominates.  Similar  conclusions  hold  v:hen  VJa-FCp  ,  r ;  A  ) . 
Case  (i)  is  of  obvious  interest.  Case  (ii)  arises  in 

estimating  the  Improvement  under  SEL  of  the  James- 
Steln  estimator  of  the  multivariate  normal  mean,  see, 
e.g.,  Efron  and  Morris  (1976). 

Two  directions  will  be  pursued.  In  the  first  a 
simple  approach  for  uniformly  improving  upon  MVU  esti 
mators  Is  described  and  illustrated.  In  the  second 
Bayesian  procedures  are  characterized  and  illustrated 
This  effort  extends  earlier  work  of  Perlman  and 
Rasmussen  (1975),  Neff  and  Strawderman  (1976),  and 
DeWaal  (197*0- 


2.  DISTRIBUTION  THEORY  AND  NOTATION 


In  the  NC  x2  case,  the  joint  density  f(Z,L,A) 
arises  from  f(Z|L,A)  =  f(ZjL)  =  x2(p+2L),  f(L(A) 

=  Po(A)  (Poisson)  and  a  prior  x(A)  on  A.  If 
A  %  y/2  x2(P)  (arising  from  A  =  uTy/2  with  y%Np(0,Yl) 
as  discussed  in,  e.g.,  Perlman  and  Rasmussen  (1975) )> 
f ( A | L, Z )  =  f ( A | L)  =  p/2  x2(P+2L),  where  p  =  y(Y+1)-1, 
f ( Y | Z )  =  p/2  x2(P>  pZ/2),  f (L)  *  NB(p/2 ,  p)  (negative 
binomial  with  mean  pp(l-p)  ^/2),  f(Z)  =  (y+1)x2(p) 
and  f ( L | Z )  =  Po(pZ/2). 

If  W  =  Z/U  where  Z  'v  x2(P">*)  independent  of 
U  %  x2(r)>  w  has  a  nonnormalized  NC  F  distribution. 

The  joint  density  f(W,L,A)  arises  from  f(W|L,A) 

=  f ( W | L )  =  r-1(p+2L)F(p+2L,r) ,  f ( L | A )  and  t(A)  as 
above . 

The  fact  that  f(Z|L,A)  -  f(Z|L),  f(W[L,A) 

=  f(W|L)  is  useful  in  improving  upon  MVU:s.  The  fact 
that  regardless  of  x,  f ( A | L , Z )  =  f ( A | L) ,  f(A|L,W) 

=  f(A|L)  is  useful  for  finding  Bayes  estimators.  In 
this  vein,  if  g(A)  is  to  be  estimated,  let  b(L) 

=  E(g(A)|L),  a(L)  be  such  that  g(A)  =  E(a(L)|A).  Then 

a(L)  helps  to  find  Improved  estimators,  b(L)  helps  to 
find  Bayes  estimators. 

In  Sections  3  and  U  the  methodology  is  examined 
through  a  collection  of  examples.  A  summary  table  is 
given  at  the  end  to  unify  the  findings. 

3-  FINDING  IMPROVED  ESTIMATORS 

We  describe  a  simple  approach  which  has  been  suc¬ 
cessful  in  creating  classes  of  estimators  that  uni¬ 
formly  improve  upon  the  MVU  of  g(A)  under  SEL.  We  then 
illustrate  with  several  examples.  For  the  NC  x2  case, 
let  T(Z)  be  the  MVU  of  g(A).  When  will  T(Z)  +  <J>(Z) 


uniformly  improve  upon  T(Z)?  The  conditional  estima¬ 
tion  problem,  estimating  a(L),  provides  sufficient  con¬ 
ditions  while  enabling  us  to  work  with  the  simpler 
central  x2  distribution.  We  note  that  T  is  unbiased 
for  a(L).  By  direct  calculation  we  can  show  that  the 
difference  in  SEL  between  T  and  T+<f>  is: 

-E(  I  (  L)  |  A  )  -  2  cov(a(L),<j>(Z)  |  A)  (1) 

where 

I ,  ( L)  =  E(4>2]L)  +  2E(  ( T— a )  <t>  |  L)  (2) 

<P 

and  the  covariance  is  with  respect  to  the  joint  distri¬ 
bution  of  L  and  Z  given  X. 

Hence  if 


I . ( L)  <  0  VL  (3a) 

4>  _ 

with  strict  inequality  for  some  L  and 

cov(a,4>)  <  0,  (3b) 

then  T+ <j)  will  dominate  T.  If  a(L)  is  monotone,  restric¬ 
tion  of  <{)  such  that  E ( <t>  |  L )  is  nonincreasing  in  a(L) 
will  satisfy  the  covariance  condition.  If  the  range 
of  g(A)  is  a  subinterval  of  r\  T+4>  should  be  similarly 

restricted.  In  particular,  if  g(A)  >_  0,  the  positive 
part  estimator,  [T  +  4>]+,  dominates  T+4>. 

For  the  NC  F  all  the  above  remarks  apply  with  Z 
replaced  by  W. 


Examples : 


<J> 


(2) 


( 1 )  8 

Two  convenient  choices  of  <J>  are  4>v  „  =  aZ  ,  «  _ 

07  (X,P  &,P 

=  aep  in  which  case  the  fact  that,  n  <  1/2,  m  >  -p/2, 


E(ZmenZ 


>m 


T  >,  _  r(p/2  +  L  +  m) _ 

r(p/2  +  L)  (l-2n)p/2  +  L  + 


m 


will  be  helpful.  We  assume  p  >  5. 


! 


\ 


1 


s 


i)  If  g(X)  =  X,  T(Z)  =  (Z-p)/2  and  b(L)  =  L.  For 
4) ^ ^  (3b)  will  be  satisfied  if  a  <  0,  M  0  or  a  >  0, 

-p/2  <  8  <  0.  For  <p^2^  (3b)  will  be  satisfied  if  a  <  0, 

0  <  6  <  1/2  or  a  >  0,  8  <  0.  Clearly,  for  any  a, 8 

cov(L,  <j>  ^  ^ )  -*•  0  as  X  -*■  0,  i  =  1,2.  (4) 

If  6  >  -p/4,  !<})(]_)(  L)  =  2Ba(T(p/2  +  L))_1 

•(2SaT(p/2  +  L  +  28)  +  2BT(p/2  +  L  +  8)).  Hence 

-p/4  <  8  <  0  requires  0  <  a  <  -2^  B8T(p/2  +  L  +  8) 

. (f(p/2  +  L  +  28) )_1  VL.  Since  the  right-hand  side  is 

smallest  at  L  =  0,  we  obtain  the  condition  0  <  a  < 

-21-B8r(p/2  +  8)(f(p/2  +  28))_1  in  agreement  with  Neff 

and  Strawderman  (p.  66)  and  including  Perlman  and 

Rasmussen  (p.  464).  If  a  exceeds  this  condition 

I  (i)(°)  >  i.e.,  with  (4),  T  +  <j>^  can't  dominate  T. 

4*  l  —  R 

If  B  >  0,  we  require  0  >  a  >  -2  P8T(p/2  +  L  +  8) 

.  ( r (p/2  +  L  +  28))_1  VL  which  is  impossible  as  L  +  ■». 

Pasmuss'-  (1973)  shows  that  when  8  =  1  no  <J>  of  the  form 

aZ  +  <5  yields  T  +  $  which  dominates  T,  i.e.,  no  linear 

estimator  can  dominate  T. 

If  6  <  1/4,  I  (2)(L)  =  a(1‘  46)"(p/2  +L) 

- 1  ^ 

•  a+28 ( 1-28 )  (p+2L)n1(L)  where  n^lL) 

=  [(1-4b)/(1-28)]P/2  +L.  Similar  to  the  4>(1)  case,  8<0 
requires  0  <  a  <  -2p8(l-28)  ^  n-^(O)  to  have  I  <  0  VL. 

If  a  exceeds  this  condition,  with  (4),  I  ^^(0)  >  0,  i.e. 

( P )  ^ 

T  +  4>v  ;  can't  dominate  T.  If  0  <  8  <  1/4,  no  a  <  0 

works  if  L  is  sufficiently  large. 

ii)  If  g(X)  =  E(Z  ''■)  (up  to  a  constant  g(X)  is  the 
improvement  of  the  James-Stein  estimator  of  the  multi¬ 
variate  normal  mean  when  the  variance  is  known),  T(Z) 

=  Z"1,  b(L)  =  (p-2  +  2L)_1.  For  4>(1)  (3b)  holds  if  a  <  0, 

-p/2  <  8  <  0  or  a  >  0,  8  >  0.  For  (3b)  holds  if 

a  <  0,  8  <  0  or  a  >  0,  0  <  6  <  1/2. 


4 


If  8  >  -p/4,  I  (1)(L)  =  2B-s(r(p/2  +L))"1 

a  $  _  i 

• ( 2paT (p/2  +L+20 )  -  20(p-2+2L)  T(p/2  +L+B-1)).  When 

p  <  8  with  B  >  -p/4,  (i.e.,  B  >  -2),  21_BB(p-2+2L)-1 
•T(p/2  +L+B-1)  (T(p/2  +L+2B))-1  >  0  as  L  >  no  T+<t>(1) 
which  dominate  T  are  revealed.  If  p  >  8,  -p/4  £  8  1  -2 
0  >  a  >  21~6B(p-2)"1r(p/2  +B-l)(r(p/2  +2B))"1,  T+<J>(1) 
will  dominate  T. 

If  B  <  1/4 ,  I  (2)(L)  =  a(l-4B)“(p/2  +L) 

i 

• (a-4B(p-2+2L)  n-,(L)).  If  B  <  -(p-4)  and  0  >  a  > 

_  i  -*■  _  1 

4B(p-2)  n1 (0)  (B  <  -(p-4)  is  needed  to  insure  that 

4b(p-2  +  2L)-1  p..(L)  is  largest  at  L  =  0),  T+4>  ^ 2  ^  will 
dominate  T.  If  0  <  B  <  1/4,  no  a  >  0  works  if  L  is 
sufficiently  large. 

iii)  If  g ( X )  =  cX  <=>  g( X )  =  ecX,  c  >  -1,  T ( Z ) 

=  (l-d)p/2edZ/2  where  d  =  c(c+l)_1,  a(L)  =  (c+l)L.  If, 
in  fact,  c  >  0  we  can  show  that  for  both  4> ^ ^  and  <J^2^ 

there  exist  a, 6  such  that  T+(j/2^  dominate  T. 

We  omit  the  details. 

iv)  If  g( X )  =  X2,  T(Z)  =  Z2/4  -  (p+2 ) Z/2 

+  p(p+2)/4  and  a(L)  =  L(L-l).  Since  a  is  monotone  on 

the  support  of  L,  (3b)  will  be  satisfied  for  and 

( 2 ) 

<p  over  the  same  ranges  as  in  example  (i). 

If  B  >  -p/4,  I  (1)(L)  =  28ct(r(p/2  +  L))_1 
a  41 

• ( 2  aT (p/2  +L+2B)  +  2 B ( B  +  2L-1 ) T (p/2  +L+B)).  It  is 

apparent  that  if  -p/4  <  B  <  0,  a  >  0,  I  ^-^(0), 

<P 

I  d)(D  >  0;  if  0  <  B  <  1>  a  <  0,  I  ^^(0)  >  Here 

$  ci)  $ 

T+$v  can't  dominate  T.  When  g  >  1,  no  a  works  if  L 
is  sufficiently  large;  this  approach  doesn't  reveal  any 
T  +  4/^  which  dominate  T. 

If  B  <  1/4 ,  I  (2)<l)  =  a(l-4B)_(p/2  +L) 

• (a+2B(l-2B)“2n, (L) (4L2(1-B)  +  2L(p+2fi)  +  Bp(p+2))). 


Again  we  do  not  reveal  any  T  +  $ 


(2) 


which  improve  upon  T. 


-2  BZ 

However  <p  =  aZ  ep  ,  for  example,  can  be  used  success¬ 
fully.  We  omit  the  details. 

Examples:  NC  F 

We  illustrate  for  <p  Q  =  aW^.  (The  more  general 

81  8p  a,B 

<t>  =  aW  (1+W)  "  can  be  used  to  broaden  the  conclusions. 
We  omit  the  details.)  It  is  useful  to  note  that, 

-p/2  <  B  <  r/2. 


e(wB|L)  *  f(p/2  +L+B)  f ( r/2  -B) 
r(p/2  +L)  r (r/2 ) 

which  increases  in  L  if  B  >  0,  decreases  in  L  if  3  <  0. 
We  assume  p,r  5- 

v)  If  g(  A  )  =  A,  T  ( V/ )  =  ( ( r-2  )W-p  )/2  and  a(L)  =  L. 
Condition  (3b)  will  be  satisfied  if  a  >  0,  -p/2  <  3  <  0 

or  a  <  0,  0  <  3  <  r/2  and  for  any  a, 3  clearly  cov(L,<}>) 

-►  0  as  A  -*■  0. 

For  -p/4  <  3  <  r/4 ,  I.(L)  =  an0(L) 

• (a+23n^(L) (p+r+2L-2) (r-23-2)  )  where 

_  /T%  _  r(p/2  +L+23)  r (r/2  -23) 

"  r(p/2  +L)  r ( r/2 ) 

n  (L)  =  r(p/2  +L+B )  f(r/2  -B) 

3  r(p/2  +L+23 )  r ( r/2  -  23) 

Note  that  no(L)  =  0(L-®)  and  that  (p+r+2L-2)n:}(L) 
increases  in  L  for  -p/it  <  3  <  p(p+r-4)  .  Hence  — p/i| 

<3<0,  0<a<  -23^(0)  (P  +  r-2)  (r-23-2)  1  or  0  <  3 
<  p(p+r-4)  *,  0  >  a  >  -23113(0)  (p+r-2)  (r-23-2)'1  yields 
T+<{»  which  dominates  T.  For  3  <  0,  if  a  exceeds  the  con¬ 
dition,  1.(0)  >  0.  The  case  3  =  -1  is  discussed  in 
Perlman  and  Rasmussen  (p.  467).  At  3  =  1,  n^(L) 

• (p+r+2L-2) (r-23-2)'1  =  (p+r+2L-2) (p+2L+2)_1  >  1  so  that 
-2  <  a  <  0,  T+<j)  will  dominate  T.  In  fact,  if  $  =  aW+y 


c 

i 

« 

« 

x 


< 


< 


a  more  general  dominating  family  of  linear  estimators 
can  be  created  in  agreement  with  Perlman  and  Rasmussen 
(p.  465-6). 

vi)  For  g(A)  =  E(W~^)  (up  to  a  constant,  g(A)  is 
the  improvement  of  the  James-Stein  estimator  when  an 
independent  estimator  of  the  variance  having  r  degrees 
of  freedom  is  used),  T  =  W  ^  and  a(L)  =  r(p-2+2L)-"'". 
Here  a  >  0,  0  <  B  <  r/2  or  a  <  0,  -p/2  <  B  <  0  will 
satisfy  the  covariance  condition. 

For  -p/4  <  B  <  r/4 ,  I  (L)  =  an2(L)  (a-4Bn3(L)ni|  (L) ) 
where  n^(L)  =  (p+r+2L-2 ) (p+2L-2 )-1 ( p+2L+2B-2 )-l  and 
n^(L)*n1J(L)  increases  in  L  if  -p/4  <  B  <  -(p(p-2  +  2r) 

• (p (p-2 ) +r (p+2 ) )-1 .  For  B  in  this  range,  if  0  >  a 
>  4Br)2 ( 0 ) n ij  ( 0 ) ,  T+(J>  will  dominate  T.  At  B  =  -1, 

n-j(L)ni,(L)  >  (r+2)  1  which  requires  0  >  a  >  -4(r+2)_1. 

•i  4  _1  _l  _i 

Hence  cW  will  dominate  W  when  (r-2)(r+2)  <  c  <  1. 


4.  BAYES  ESTIMATORS 


In  developing  Bayesian  procedures 
again  turn  to  the  conditional  problem, 
relation 


.L  -A 

tt(L)  =  —  t  (  A  )dA 

L! 


under  SEL, 

i  .  e  .  ,  Z/L . 


we 

The 


(5) 


shows  which  tt(L)  can  arise  as  priors.  In  fact,  since 
(5)  is  an  instance  of  the  classical  moment  problem,  if 
-rr(L)L.'  is  a  "moment  sequence"  (see,  e.g..  Feller  (1966) 
Sec.  VII.  3  Tor  conditions)  tt(L)  uniquely  determines 
t(A).  A  useful  case  is,  v  >  -  , 


A  %  G(p/2  +v;y  ^)  (Gamma  with  mean  (p/2  +v)y),  (6) 

L  NB(p/2  +v;  p)  (p  =  Y(y  +  1)"1)  . 

Under  (6),  A|L  'v  G(p/2  +v+L;  p_1). 


4 


7 


I 

i 

P 


E 


Recalling  that  b(L)  =  E(g(A)|L),  for  the  NC  x2 
we  have  by  direct  calculation  E(g(A)|Z)  =  E(b(L)|Z), 
i.e.,  we  can  calculate  Bayes  rules  using  the  central  x2 
distribution.  The  same  conclusions  hold  for  the  NC  F 
with  W  replacing  Z. 

Examples:  NC  x2 

For  a  particular  tt  the  Bayes  rule 

<5  (Z)  =  E  ( b  ( L )  |  Z )  =  -b  ill  PL2  +L) ) >  ( 7 

v  Z(Z/2)Ltt(L)  (  T (p/2  +L)  )-i 

Denote  the  denominator  in  (7)  by  J  (Z). 

TT 

i)  If  g(A)  =  A,  then  from  (5)  b(L)  =  (L+l) 

•  tt  (  L+l )  (  tt  (  L )  )~ 1 . 

Using  this,  straightforward  manipulation  yields 


6  (Z)  =  ?Z  ■--- - +  p  — - 

"  J  (Z)  J  (Z) 

TT  TT 


Setting  6  (Z)  to  be  the  MVU  of  X  in  (8)  leads 

TT 

second  order  homogeneous  linear  differential 
whose  general  solution  is 


(8) 


to  a 

equation 


J^Z)  =  e  Z/2(c1  +  c2/z  p/2  eZdZ)  . 


(9) 


By  the  definition  of  J  c „  cannot  equal  0.  But  then 

1T  ^  z/2 

after  multiplying  both  sides  of  (9)  by  e  and  dif¬ 
ferentiating,  we  see  that c^  ?  0  can't  work  either. 

Thus  the  MVU  can't  be  the  limit  of  Bayes  or  extended 
Bayes  clarifying  its  inadmissibility. 

Applying  the  NB  priors  in  (6)  to  (8)  and  denoting 

the  resulting  rule  by  6  (Z),  we  have 

v ,  p 


6 

v ,  p 


(Z) 


p  (  2"  +  v  + 


ZJ 

IT 

'  J 

TT 


)  • 


(10) 


8 


It  is  straightforward  to  show  that  for  fixed  p. 


ZJ 


<  6 


.if  v,  < 
Vo  P  1 

> 


At 


6n  =  p(pZ+p)/2  which  is  dis- 
u  j  p 


increases  in  v  whence  6, 

v  =  0,  J^Z)  =  cepZ/2, 
cussed  at  length  (particularly  when  p  =  1)  in  Perlman 
and  Rasmussen.  Corresponding  to  the  noninformat ive 
prior,  tt  ( L )  =  1  from  t(A)  =  1  (which  arises  in  (6)  at 
v  =  -  (p-2  )/2 ,  p  =  1),  we  have  &  ? 

2 

+  2  so  that  the  MVU,  T  <  6  n  Vv. 

V  J  -L 


( j  ( z ) r(p/2  -D) 


-i 


In  fact ,  6 


-(^-),1 


=  T+4»  with  4>  satisfying  (3b),  but  we 


are  unable  to  show  that  this  6  dominates  T.  We  note 
that  Neff  and  Strawderman  derive  a  sublcass  of  (8)  of 
proper  Bayes  estimators  for  A  arising  from  a  twc-stage 
prior  distribution  and  show  that  none  of  these  dominate 
T. 

ii)  Generally  from  (5)  if  g(A) 


integer  b(L) 

2 

for  A  the  Bayes  rule  becomes 


(L+r)r  it (L+r )  (tt(L)  ) 


-1 


A  ,  a  positive 
In  particular. 


jf'v)(z) 


( '* » i 

;(z) 


J(,,)(z) 


VZ>=/4Z2  ~^( Yj~+liZ('P+2')  j^(z)  +p^p+?)  j^(z)  ’ 

Using  an  argument  similar  to  that  below  (9),  it  may  be 
possible  to  show  that  the  MVU  can't  be  the  limit  of 
Bayes  or  extended  Bayes. 

to  (11),  we  obtain,  for  example,  v  =  0,  p 
+  p3  Z(p+2)/2  +  p2  p(p+2)/H. 

iii )  If  g( A )  =  ecA,  under  (6)  with  c  <  p  . 
b(L)  =  (l-pc)_  p,/2  +V+L^.  Bayes  rules  take  the  form 


Applying  the  NB  priors  in  (6) 
h  O 


-1 


6„(Z) 


•  (l-pc)"(p/2  +v)J  ((1-pc)  hl/J  (Z>  .  (12) 


9 


At  v  =  0  we  obtain  (1-pc)  p//2  ep  c Z / 2 ( 1  pc).  The  ^na<j_ 
missible  MVU  estimator  would  arise  if  p  could  equal  -1. 

iv)  If  g  (  X )  =  ECZ'1)  =  E( (p-2+2L)-1 | X) ,  then 
under  (6),  b(L)  =  E( (p-2+2M)-1 | L)  where  M|L  ^  NB 
(p/2  +v+L),  p(p+l)  1),  i.e.. 


,  /T  x  _  v  r(p/2  +v+L+m)  f  p 

b(L)  mlQ  f (p/2  +v+L)mf  ^p~2  2  ^  ^p+1 


.m,  1  ^p/2  +v+L 
J  lp+l; 


When  v  is  a  nonnegative  integer,  we  have  the  identity 


*■  r(p/2  +v+L+m)  ,  n  ,-lf  p  xp/2  +m-l 

2  r ( p / 2  +v+LTmT  (p-2+2m)  (p+l} 


I  (^llP-Wnr1  Pp/2  tn-1 
n  =  0 


(derivable  by  considering  the  indefinite  integral  with 
respect  to  p  of  pP//2  -2(p+l)L+V  directly  and  through 

Its  equivalent  negative  binomial  expansion) .  This 
yields  the  Bayes  rule  in  the  form 


(l-p)p/2(p+l)~(v*I}  .  “  f(p/2)  +  L+v)  (  p  ZvL 

^p+1  2 ) 


J  (Z)f(p/2  +v) 


L=0  r(p/2  +L)L! 


(13) 


v+L  , 

•  l  (V^L)(p-2+2n)~1  pn 
n=0  n 


At  v  =0,  interchanging  order  of  summation  in  (13),  we 

obtain  the  rule  as  (p+1)  bE((p-2+2L)  b|Z)  =  (p+1)  b 
2 

•g(p  Z/2(p+l))  since  L|Z  has  a  Poisson  distribution. 


Examples:  NC  F 

For  g(X),  b(L),  under  it,  the  Bayes  rule  becomes 
(in  terms  of  V  =  W(l+W)  ^  ^  NC  Beta) 


6ii(V)=E(a(L)  |  V)  = 


lb  (  L)  VLti  (  l)  r(p^-  +L)(r(|  +L) )  1 
ZVL7r(L)f(E^  +L)(f(|  +L))"1 


( 1^  ) 


Is 


« 


X 


B 


i 


4 


I 


Denote  the  denominator  in  (lA)  by  K  (V). 

71 

v)  If  g(A)  =  A  it  is  straightforward  to  show  that 


_(E±I1)  P+r 

67r(V)  =  (K7i(V))~1V  2  /V  2 


-1 


(VKrWfK^m).  (i5) 


Under  (6),  (15)  becomes  (denoting  the  rule  by  6  ) 

v,p 


6'..  n(V)  =  p[p/2  +v+VK  /K  ]. 

V  >  P  TT  TT 


(16) 


Note  the  similarity  between  (16)  and  (10).  It  is 
straightforward  to  show  that  for  fixed  p.  v  +  VK  /K 

TT  IT 

increases  in  v  whence  6  <6  if  v-,  <  v~.  At 

v,  p  v~  p  1  v2 

±  y  C.  y 

V  =  o,  Krr(v)  =  (i-p)p/2(r(p/2) )-1r(£|^) (i-Pv)_(p+r)/2 

so  that  6q  (V)  =  p/2 (  (p+r )  pV(  1-pV)  '*‘  +  p)  which  is  dis¬ 
cussed  in  Perlman  and  Rasmussen  (p.  ^66),  particularly 
when  p  =  1,  (6q  ^  3  ( (p+r ) W+p )/2 ) .  Corresponding  to  the 

noninforinative  prior,  we  obtain  6  _  =  (r-2)W/<i 

-(P~),l 

p/2  +  ( W+l )  •  (  2  +  r  ( p/2~  - 1 )  •K7i(W(  1 +WT-T  )  whence 

T  <  6  i ,  Vv,  where  T  is  the  MVU  estimator. 

v,  1  ’ 

vi )  If  g(A)  =  E((p-2  +  2L)  ^  |  A  ) ,  under  (6)  calcula¬ 
tions  analogous  to  those  leading  to  (13)  yield  at  v  =  0 
the  rule  (p+1)  ^E((p-2+2L)  ^|W)  as  in  the  NC  x2  case 
except  L|W  has  a  negative  binomial  distribution. 


5.  SUMMARY  AND  DIRECTION  FOR  FUTURE  WORK 

We  present  a  summary  table  of  the  disparate  findings 
in  this  article. 

Two  problems  requiring  further  investigation  are 
the  following.  First  it  would  be  useful  to  link  the  two 
estimation  approaches  discussed  herein.  More  precisely, 
writing  the  Bayes  estimators  in  Section  4  in  the  form 


< 


11 


T+4>,  where  T  is  the  corresponding  MVU,  can  we  use  the 
method  of  Section  3  (or  some  other  argument)  to  obtain 
Bayes  rules  which  uniformly  improve  upon  the  MVU?  Do 
any  Bayes  rules  dominate  the  MVU?  Second,  under  the 
priors  in  (6)  at,  for  example,  v  =  0  the  marginal  dis¬ 
tribution  of  Z  is  (l-p)_1X2(p)  suggesting  that  for  the 
NC  x2  case  convenient  "empirical  Bayes"  estimators  can 
be  developed.  Do  any  of  these  estimators  dominate  the 
corresponding  MVU?  Can  the  method  of  Section  3  help  in 
this  regard? 

Finally,  the  meticulous  and  provocative  effort  of 
the  referee  is  acknowledged.  The  paper  is  much  improved 


TABLE  I 


gU) 


Improving 
Estimator  of  the 
Form  T+4> 


Bayes  Estimator 
under  SEL 


~  aZ 


B 

BZ 


4>  =  ae 
See  Ex.  3(i) 


General  form:  (8) 

Under  priors  in  (6): (10) 
See  Ex.  4(i) 


*  =  aZ2egZ 
See  Ex.  3(iv) 


General  form:  (11) 
See  Ex.  *J(ii) 


or  e 


CA 


<p  =  aZ 


e 

BZ 


<})  =  ae 
See  Ex.  3(iii) 


Under  priors  in  (6): (12) 


See  Ex.  H(iii) 


E(Z_1) 


4>  =  aZ6 

.  6Z 

4>  =  ae 

See  Ex.  3(ii) 


Under  priors  in  (6):  (13) 


See  Ex.  4(iv) 


B 


4>  =  aW 


See  Ex.  3(v) 


General  form:  (15) 

Under  priors  in  (6): (16) 
See  Ex.  4(v) 


E ( W~ 1 ) 


<p  =  aWB 
See  Ex.  3(vi) 


See  Ex.  4(vi) 
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